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TECHNICAL MEMORANDUM 


GEOMETRIC PROGRAMMING PREDICTION OF DESIGN TRENDS 
FOR OMV PROTECTIVE STRUCTURES 

I. INTRODUCTION 


The orbital maneuvering vehicle (OMV) design activity presents many challenges to the 
design engineer. The OMV will be used to transfer Earth-orbiting satellites to different orbits and 
will then be returned to Earth aboard the space shuttle orbiter. The design ot an optimized protec- 
tive shield which can prevent orbital debris particles from impacting and damaging the OMV could 
be a major design and development effort. Thousands of large orbiting man-made particles are 
being tracked each year by the space science community; many more smaller untrackable particles 
are known to exist from evidence of impact damage on returned satellites. The OMV, with 
exposed high-pressure fuel bottles, must be adequately protected for its average orbital lifetime of 7 
days to prevent possible penetration and catastrophic rupture of the bottles. 

As for any launched space vehicle, weight is a primary structural concern. The geometric 
programming optimization technique has been employed to optimize the metallic honeycomb struc- 
ture, which is currently required as a stiffening panel on the OMV configuration. I his technique 
determines the facing sheet and core areal densities needed to meet a protection requirement at the 
minimum structural weight. 

This report covers the first step in using geometric programming to optimize the honeycomb 
panel configuration. By parametrically relating the material design variables in model forms 
indicated from past experience, general trends can be predicted for the honeycomb panel as an 
optimized protective structure. This information can then be used as a stepping stone for further 
work in this area by providing researchers with clearer insight into the relationships of the design 
parameters. 


II. PROBLEM FORMULATION 


The formulation of the protective structures design optimization problem begins with the 
specification of an objective function in the form of engineering models and the statement of 
applicable constraints in terms of the design variables. For the minimization of protective structures 
weight, it has been shown that an objective function in the form of structural mass per unit area is 
sufficiently representative, provided the spacecraft configuration to be protected is not dominated by 
high curvature surfaces [1,2,3]. This requirement is met by OMV. For the honeycomb configura- 
tion shown in figure 1 , the system mass per unit area is given by 
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Figure I. Notation for honeycomb structure configuration. 


The wall thickness required to resist penetration/spallation is a function of the impact parameters, 
structural material properties, and structural configuration f 1 ,31 - For the purposes of this develop- 
ment, we assume a single term posynomial (or monomial) relationship given by 


h = K^ P y:dY 6 v\ M (Q))y 2 '' 


( 2 ) 


We will see later that this assumption is not particularly restrictive and can be greatly generalized 
with equivalent results. Five different problem formulations are considered by substituting equation 
(2) into equation (1). These formulations involve different combinations of independent variables. 
Case (i) considers 7, and t c as independent variables. Case (ii) also includes the material densities 
of the facing sheets and core. Case (iii) is the same as case (ii) but without the second facing sheet 
material density as an independent variable. Case (iv) is the same as case (iii) with some dimen- 
sional analysis considerations to reduce the number of required estimated parameters. Finally, case 
(v) is the same as case (iv) with the second facing sheet material density independent. 

Iii. STRUCTURAL OPTIMIZATION USING GEOMETRIC PROGRAMMING 


The five cases are solved by maximizing the corresponding dual problem subject to the 
derived linear constraints. The dual objective function takes the product form as specified by the 
arithmetic-geometric inequality [4-7], The linear constraints include normality as well as orthog- 
onality constraints. The normality constraint merely specifies that the dual variables sum to 1. This 
is a required condition for the application of the arithmetic-geometric inequality. The orthogonality 
conditions specify for each independent variable that the product of the appropriate dual variable 
with the independent variable’s exponent summed over all of the terms in the primal posynomial be 
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equal to 0. The main reason for transforming the problem to a dual form is to take advantage of 
the linear constraints in a low degree of dilficulty problem. Additionally, the dual problem results 
in analytic solutions for zero degree of difficulty problems as given here. Finally, geometric 
programming provides globally optimal solutions for problems in posynomial form [4—7], For case 
(i), the dual problem is given by 
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The solution is given by 
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Since the dual variables must all be positive, we have 


a 2 < 0 a 4 < 0 . 


(7) 


Objective function solution and conversion to the primal variables will be deferred until case (v). 
For case (ii), the linear constraints are given by 
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This system is inconsistent unless 
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If condition (10) did not hold, then each of the dual variables would be identically zero. This 
would contradict the normality condition that the dual variables must sum to 1 If condition (10) 
holds then 


ci\ = a 2 a 3 = a 4 
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and (6) again provides the solution for the dual variables. 

Case (iii) results in the same solution as case (ii) without requirement (10). In case (iv), the 
structural mass per unit areas are normalized by the projectile diameter density product to assure 
dimensionless variables and reduce the amount of regression required. The functional relationship is 
given by 
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with the dual problem given by 
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The dual variable solution is again given by (6). Case (v) includes second facing sheet density as 
an independent variable. The results are identical to case (iv) with requirement (10) included. The 
minimum system mass per unit area for the structure is found by substituting (6) into (13) and is 
given by 


W 0 = 


“2 


( a 2 + a 4 - 0 

* 1 

f K(a 2 + a 4 - 1) 

^ J 


l "I 


-1 



\ a 2+ a 4~ ] 

) 


(16) 


where 
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Additionally, the optimal mass per unit areas for the first facing sheet, the second facing sheet, and 
the honeycomb core are found by multiplying the minimum system mass per unit area by the 
corresponding dual variable. These values are given by 
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IV. PARAMETRIC DESIGN TRENDS 


This section gives parametric trends for solutions (16) to (20) in terms of the impact system 
constant, a \ ( = a 2 ) and a 3 ( = a 4 ). Figure 2 shows the optimal mass per unit areas as functions of 
the impact system constant for a 2 = -0.5 and a 4 = -0.5. The First sheet and honeycomb core 
lines are equal and bounded above by the second sheet mass per unit area. 

Figure 3 shows the decreasing relationship between minimum system mass per unit area and 
a 4 , the honeycomb core mass per unit area exponent. Several curves are shown, corresponding to 
different values of a\ ( = 02 )- The case a f =0 corresponds to a hypothetical structural configuration 
without a first sheet material. As the honeycomb core exponent increases, the core is less effective 
in terms of penetration resistance, and the structural configuration increasingly becomes a bumper/ 
wall situation. The crossings of the a, = -0.75 and a, = -1.0 lines represent a relative tradeoff 
in the utility of the honeycomb core and the penalty imposed by its weight. 

Figure 4 (fig. 6) shows the nonmonotonic relationship between optimal first (second) sheet 
mass per unit area and a 4 for various values of a 2 . Increases in a 4 imply a decreasing role in the 
honeycomb core, which requires increases in the first and second sheet mass per unit areas, to a 
point. Beyond that point (corresponding to a 4 = —0.25), the honeycomb structure becomes increas- 
ingly inefficient relative to a bumper/wall configuration. At this point, decreases in first and second 
sheet mass per unit areas are in order. 


Figure 5 shows the sharp decreases in optimal honeycomb core mass per unit area as a 
function of a 4 . Note that the anomalous curve a x — 0 is bounded by a\ = -0.75 and a\ = -1 .0. 

Figures 7, 8, and 9 show the relationships of optimal ratios for the first, second, and 
honeycomb core materials as a function of a 4 for various values of a \ . This ratio is defined as the 
optimal mass per unit area of the material divided by the minimum system mass per unit area and 
is given by the dual variable values given by equation (6). Note that for the first and second 
sheets, this ratio is monotonically increasing, while for the honeycomb core, the ratio is monotoni- 
cally decreasing. 
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Minimum System Mass Per Unit Area 


First Sheet Mass Per Unit Area Exponent = -0.5 
Honeycomb Core Mass Per Unit Area Exponent = -0.5 



System 
First Sheet 
Honeycomb Core 
Second Sheet 


Figure 2. Optimal mass per unit areas versus impact system constant ( K ). 
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Figure 3. Minimum system mass per unit area versus honeycomb core mass 

per unit area exponent ( 04 ). 





Optimal Honeycomb Core Mass Per Unit Area 
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Figure 4. Optimal first sheet mass per unit area versus honeycomb core mass 

per unit area exponent ( a 4 ). 




Impact System Constant = 1 .0 



Figure 5. Optimal honeycomb core mass per unit area versus honeycomb core mass 

per unit area exponent ( a 4 ). 
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Optimal First Sheet Ratio n Optimal Second Sheet Mass Per Unit Area 
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Optimal second sheet mass per unit area versus honeycomb core mass 
per unit area exponent (u 4 ). 
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Figure 7. Optimal first sheet ratio versus honeycomb core mass 
per unit area exponent ( a 4 ). 
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Optimal Second Sheet Ratio Optimal Honeycomb Core Sheet Ratio 
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Figure 8. Optimal honeycomb core ratio versus honeycomb mass 
per unit area exponent (a 4 ). 
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Figure 9. Optimal second sheet ratio versus honeycomb core mass 
per unit area exponent (tf 4 ). 
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V. CONCLUSIONS 


Using the geometric programming optimization technique, the following conclusions can be 
drawn for a honeycomb-type panel: 

1 . Optimal ratios are monotonically increasing for the first and second facing sheets and 
monotonically decreasing for the honeycomb core. 

2. Optimal mass per unit areas are unimodal for the first and second facing sheets and 
monotonically decreasing for the honeycomb core. 

3. The numerator of the impact system constant is not restricted to equation (17). Any form 
for the numerator will result in identical results for the parametric trends given in figures 1 to 8 
and equations (16), (18), (19), and (20). 

4. The minimum system mass per unit area is an essentially monotonically decreasing func- 
tion of honeycomb core mass per unit area exponent. 

5. The optimal mass per unit area of the facing sheets, honeycomb core, and system are 
monotonically increasing functions of the impact system constant. 

In addition to these conclusions, one more point should be made. Geometrical con- 
siderations, particularly involving honeycomb cell size and shape, have been purposely omitted to 
reduce the number of estimated parameters required for model development. Certainly, equation 
(17) can be generalized to include these parameters, provided the impact test database supports this 
development. An interesting tradeoff involving cell size, penetration resistance, and cell density 
could result. 

The conclusions from this study indicate that the geometric programming method can aid 
the designer of a honeycomb panel for optimal protection from particle impacts while minimizing 
weight. Additionally, an important conclusion is that the geometric programming method of 
optimization can be helpful in many fields of design engineering where optimal solutions may 
involve tens or hundreds of independent variables and design parameters. These design problems 
may be solved with geometric programming more efficiently than with other methods, and result in 
globally optimal, rather than locally optimal, solutions. 


11 


REFERENCES 


1. Mog, R.A.: “Spacecraft Protective Structures Design Optimization.” 28th Aerospace Sciences 
Meeting, AIAA 90-0087, January 1990. 

2. Mog, R.A., Lovett, J.N., and Avans, S.L.: “Global Nonlinear Optimization of Spacecraft 
Protective Structures Design.” NASA TM- 100387, January 1990. 

3. Mog, R.A., and Price, D.M.: “Optimization Techniques Applied to Passive Measures for 
In-Orbit Spacecraft Survivability - Final Report.” SAIC HV410-12, Contract NAS8-37378, 
November 1987. 

4. Duffin, R.J., and Peterson, E.L.: “Geometric Programming With Signomials.” Journal of 
Optimization Theory and Applications, Vol. 11, January 1973. 

5. Gottfried, B.S., and Weisman, J.: “Introduction to Optimization Theory.” Prentice-Hall, Inc., 
Englewood Cliffs, New Jersey, 1973. 

6. Peterson, E.L.: “A Handbook of Geometric Programming.” Northwestern University, Evanston, 
Illinois, July 1974. 

7. Wilde, D.J.: “Globally Optimal Design.” John Wiley and Sons, Inc., 1978. 


12 



BIBLIOGRAPHY 


Adeli, H., and Kamal, O.: “Efficient Optimization of Space Trusses.” Computers and Structures, 
Vol. 24, No. 3, 1986. 

Avriel, M.: “Nonlinear Programming Analysis and Methods.” Prentice-Hall, Inc., Englewood 
Cliffs, New Jersey, 1976. 

Bazaraa, M.S., and Shetty, C.M.: “Nonlinear Programming.” John Wiley and Sons, Inc., 1979. 

Dembo, R.S.: “Current State of the Art of Algorithms and Computer Software for Geometric 
Programming.” Journal of Optimization Theory and Applications, Vol. 26, October 1978. 

Dembo, R.S.: “Dual to Primal Conversion in Geometric Programming.” Journal of Optimization 
Theory and Applications, Vol. 26, October 1978. 

Dembo, R.S.: “Sensitivity Analysis in Geometric Programming.” Journal of Optimization Theory 
and Applications, Vol. 37, May 1982. 

Ecker, J.G., et al.: “A Modified Reduced Gradient Method for Dual Posynomial Programming.” 
Journal of Optimization Theory and Applications, Vol. 26, No. 2, October 1978. 

Lavi, A.: “A User Manual for Geometric Programming.” Camegie-Mellon University, Pittsburgh, 
Pennsylvania, December 1974. 

Mancini, L.J., and Wilde, D.J.: “Interval Arithmetic in Unidimensional Signomial Programming.” 
Journal of Optimization Theory and Applications, Vol. 26, No. 2, October 1978. 

McCormick, G.P.: “Nonlinear Programming Theory, Algorithms, and Applications.” John Wiley 
and Sons, Inc., 1983. 

Mog, R.A.: “Comparison of Nonlinear Optimization Techniques for Spacecraft Protective Structures 
Design - White Paper.” SAIC, June 1988. 

Mog, R.A.: “Discrete Nonlinear Optimization of Spacecraft Protective Structures Design - White 
Paper.” SAIC, August 1988. 

Mog, R.A., and Price, D.M.: “Linear Programming Estimation of the Nonlinear Spacecraft Protec- 
tive Systems Design Optimization Problem - White Paper.” SAIC, March 1988. 

Morris, A.J.: “On Condensed Geometric Programming in Structural Optimization.” Computer 
Methods in Applied Mechanics and Engineering, Vol. 15, August 1978. 

Peterson, E.L.: “Fenchel’s Duality Theorem in Generalized Geometric Programming.” Northwestern 
University, Evanston, Illinois, October 1976. 


13 



Peterson, E.L.: “Optimality Conditions in Generalized Geometric Programming.” Northwestern 
University, Evanston, Illinois, May 1976. 

Ratner, M., et al. : “Solving Geometric Programs Using GRG-Results and Comparisons.” Stanford 
University, California, May 1976. 

Reklaitis, G.V.: “Investigation of the Computational Utility of Geometric Programming Formulation 
and Algorithms.” Purdue University, West Lafayette, Indiana, January 1981. 

Rijckaert, M.J., and Martens, X.M.: “Bibliographical Note on Geometric Programming.” Journal 
of Optimization Theory and Applications, Vol. 26, No. 2, October 1978. 

Rijckaert, M.J., and Martens, X.M.: “Comparison of Generalized Geometric Programming Algo- 
rithms.” Journal of Optimization Theory and Applications, Vol. 26, No. 2, October 1978. 

Rosenberg, E.: “A Globally Convergent Condensation Method for Geometric Programming.” 
Stanford University, California, October 1979. 

Sarnia, P.V.L.N., et al.: “A Comparison of Computational Strategies for Geometric Programs.” 
Journal of Optimization Theory and Applications, Vol. 26, No. 2, October 1978. 

Vanderplaats, G.N.: “Numerical Optimization Techniques for Engineering Design with Applica- 
tions.” McGraw-Hill, Inc., 1984. 

Whittle, P.: “Optimization Under Constraints.” John Wiley and Sons, Ltd., 1971. 


14 






NASA 

Hilt- * >. i( A> 'i * k )i ita , , ti k | 
' >4 1 H * •* 1m hi h Jr ilk H’ 


Report Documentation Page 


1. Report No 

NASA TM-4206 


2. Government Accession No. 


4. Title and Subtitle 


Geometric Programming Prediction of Design Trends 
for OMV Protective Structures 


7. Author(s) 

R.A. Mog* and J.R. Horn 


9. Performing Organization Name and Address 


George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


12. Sponsoring Agency Name and Address 


National Aeronautics and Space Administration 
Washington, D.C. 20546 


3. Recipient's Catalog No. 


5. Report Date 

June 1990 


6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 

M-632 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared by Structures and Dynamics Laboratory, Science and Engineering Directorate. 
*Science Applications International Corp., Huntsville, Alabama. 


16. Abstract 


The global optimization trends of protective honeycomb structural designs for spacecraft 
subject to hypervelocity meteoroid and space debris impacts are presented. This nonlinear 
problem is first formulated for weight minimization of the orbital maneuvering vehicle (OMV) 
using a generic monomial predictor. Five problem formulations are considered, each dependent 
on the selection of independent design variables. Each case is optimized by considering the 
dual geometric programming problem. The dual variables are solved for in terms of the 
generic estimated exponents of the monomial predictor. The primal variables are then solved 
for by conversion. Finally, parametric design trends are developed for ranges of the estimated 
regression parameters. Results specify nonmonotonic relationships for the optimal first and 
second sheet mass per unit areas in terms of the estimated exponents. 


17. Key Words {Suggested by Author(s)} 

geometric programming, honeycomb, meteroid, 
space debris, hypervelocity, 
nonlinear optimization 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category: 18 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21 . No. of pages 

22. Price 

Unclassified 

Unclassified 

20 

A03 


NASA FORM 1626 OCT 86 

For sale by the National Technical Information Service, Springfield, VA 22161-2171 

NASA-Langley, 1990 





